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THE LIKELY ORBITAL PERIOD OF THE ULTRACOMPACT LOW-MASS X-RAY BINARY 2S 0918-549 
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ABSTRACT 

We report the discovery of the likely orbital period of the ultracompact low-mass X-ray binary 
(LMXB) 2S 0918—549. Using time-resolved optical photometry carried out with the 8-m Gemini 
South Telescope, we obtained a 2.4-hr long, Sloan r' light curve of 2S 0918—549 and found a periodic, 
sinusoidal modulation at 17.4±0.1 min with a semiamplitude of 0.015±0.002 mag, which we identify 
as the binary period. In addition to 4U 0513—40 in the globular cluster NGC 1851 and the Galactic 
disk source 4U 1543-624, 2S 0918-549 is the third member of the ultracompact LMXBs that have 
orbital periods around 18 min. Our result verifies the suggestion of 2S 0918—549 as an ultracompact 
binary based on its X-ray and optical spectroscopic properties. Given that the donor in 2S 0918—549 
has been suggested to be either a C-0 or He white dwarf, its likely mass and radius are around 0.024- 
0.029 M and 0.03-0.032 R Q , respectively, for the former case and 0.034-0.039 M Q and 0.033-0.035 
Rq for the latter case. If the optical modulation arises from X-ray heating of the mass donor, its 
sinusoidal shape suggests that the binary has a low inclination angle, probably around 10°. 
Subject headings: binaries: close — stars: individual (2S 0918—549) — stars: low-mass — stars: 
neutron — X-rays: binaries 



1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) constitute a large ^ 
fraction of bright X-ray sources (Lx ~ 10 36 erg 
s _1 ) in the Galaxy. These binary systems consist of 
an accreting compact star, either a neutron star or 
black hole, and a Roche lobe-filling, low-mass com- 
panion. Thus far approximately ~200 LMXBs are 
known. Among them, there is a class called ul- 
tracompact binaries. Different from the majority of 
LMXBs which contain ordinary, hydrogen-rich mass 
donors, ultracompact systems are believed to consist of 
extremely low-mass, either hydrogen-poor or degener- 
ate, companion stars dNelson Rappaport. fc Joss! 119861 : 
lYungelson. Nelemans. fc van den Heuvell 120021 ) . As a 
result, while ordinary LMX Bs have a minimum or- 
bital p eriod around 80 min (iPaczynski fc Sienkiewiczl 
119811: iRappaport. Joss, fc Webbinkl I1982D . ultracom- 
pact systems can evolve to extraordinarily small 
binary separations with orbital periods as short 
as a few minutes (IPodsiadlowski. R appapo rt. fc Pfahll 
I2002t INelson fc Rappaportl 12003ft . These ultracompact 
LMXBs, along wit h their white d warf analogues (the AM 
CVn binaries; see lWarnerlll995f) . represent extreme and 
exotic endpoints in binary and stellar evolution. 

While the ultracompact systems had initially been as- 
sumed to be relatively rare, the number known has more 
than doubled to 11 (including 4 globular cluster sources) 
over the past few years, w ith a range of orbital peri- 
ods f rom 11 to 55 minutes jja fc Lil 120091; IZurek et all 
2009). It is likely that there are more such binaries, be- 
cause a few candidate systems have been identified either 
by their peculiar X-ray and/or optica l spectral features 
(jJuett, Psaltis. fc Chakrabartvl 120011: INelemans et al.l 
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120041 : Wang 20041) or thro ugh their unusually low optical 
to-X-ray flux ratios (iDeutsch. Margon. fc Anderso: 
2000; Bassa et al.ll2006l: lin't Zand. Jonker. fc Markward 
2007|) . To fully study and understand the ultracom- 
pact LMXB population, verification of those candidate 
systems are warranted. It has shown that the indi- 
rect methods for ultracompact binary identification may 
not be reliable dShahbaz et al.l [2007). In order to ver- 
ify their ultracompact nature, time-resolved photometry 
for detecting orbital periodic signals is needed. More- 
over once ultra-short orbital periods are found, proper- 
ties of the binary systems c an be further estimated (e.g. 
I Wang fc Chakrabartvl l200l . helping our understanding 
of these systems. In an effort to verify the ultracompact 
nature of the proposed candidates, we have undertaken 
optical observations aiming to detect orbital flux modu- 
lations. We have successfully fo und the orbital period 
of th e candidate 4U 1543-624 (jWang fc Chakrabartvl 
l200l . In this paper we report our discovery of the likely 
orbital period of another candidate 2S 0918—549. 

The LMXB 2S 0918-549 has been a bright X-ray 
source (Lx ~ 10 36 erg s" 1 ) and detected by all ma- 
jor X-ray satellites (jJuett fc Chakrabartvl 120031 and ref- 
erences therein) . On the basis of comparison of its X-ray 
spectrum to that of the known ultracompact LMXB 4U 
1626—67, the source has been suggested to be an ultra- 
compact binary w i th a neon-enriched degener ate donor 
(jJuett et al.l [20011: IJuetFfc Chakrabartvl l200l . Proba- 
bly because the binary has a low inclination angle (gen- 
erally i < 60°; iFrank. King, fc Raindl2002t ). no orbital 
signals were found in X-ray observations of the source 
(jJuett fc Chakrabartvl 12003). Th e optical counterpart 
to 2S 0918-549 was identified by iChevalier fc Ilovaiskvl 
(|1987ft . V = 21, B - V = 0.3. Based on its optical- 
to-X-ray flux ratio, the orbital pe riod of the binary has 
been suggested to be < 60 min (jJuett fc Chakrabartvl 
2003). The source distance is probably 4.1-5.4 kpc, esti- 
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mated from type-I X-r ay bursts detected from the source 
(jin't Zand et al.|[2005l ). 

2. OBSERVATIONS AND DATA REDUCTION 

Time-resolved imaging of 2S 0918—549 was carried out 
on 2008 December 5 using the 8-m Gemini South Tele- 
scope. The instru ment was Gemini Multi-Object Spec- 
trograph fGMOS: iHook et al.ll200"l . whose detector ar- 
ray consists of three 2048x4608 EEV CCDs. We used 
only the middle CCD chip (CCD 02) for imaging. The 
pixel scale is 0.073"/pixel, while the detector was 2x2 
binned for our observation. A Sloan r' filter with the 
central wavelength at 6300 A was used. We obtained 
179 continuous frames with an exposure time of approx- 
imately 24.5 sec. Including 24 sec readout time for each 
exposure, the total length of the observation was ap- 
proximately 2.4 hrs. The observing conditions during 
the early part of our observation (first 59 frames) were 
relatively poor, with the average seeing [FWHM of the 
point-spread function (PSF) of the images] being ~ 1.0" 
and a few frames having as large as 1.3" -1.5" seeing. 
For the remaining part, the conditions were good with 
most of the frames having < 0.8" seeing. 
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Fig. 1.— Gemini South r' image of the 2S 0918-549 field. Object 
X is the optical counterpart to 2S 0918—549. The nearby star 
labeled as N is 1.5" away from the target. The star labeled as C 
was used as a check star. 

We used the IRAF packages for data reduction. The 
images were bias subtracted and flat fielded. The bias 
and flat frames were from GMOS baseline calibrations 
taken during the same night. 

We performed PSF-fitting photometry to obtain 
brightnesses of our target and other in-field stars, with a 
photo metry program D0PH0T (Schcchtcr. Mat eo, fc Sahal 
1993) used. A finding char t of the target field is shown i n 
Figure [TJ As identified by iChevalier fc Ilovaiskvl (|1987l ). 
there is a nearby star (labeled as N) 1.5"away from our 
target. To avoid possible contamination from this nearby 
star caused by the poor observing conditions during the 
early part of our observation, we positionally calibrated 
the first 59 frames to a reference image that was com- 
bined from three high-quality frames. The positions of 
star TV and 2S 0918—549 were determined in the ref- 
erence image, and were fixed at the positions for pho- 
tometry of the first 59 frames. At last we excluded 



three frames among them from the data. From the three 
frames, the brightness measurements of star N and other 
in-field stars obtained were not consistent with the av- 
erage brightnesses from other frames, and we note that 
the three frames have the seeing of 1.3" 1.5". 

Differential photometry was performed to eliminate 
systematic flux variations in the images. Three isolated, 
nonvariable bright stars in the field were used. The 
brightnesses of our targets and other stars in each frame 
were calculated relative to the total counts of the three 
stars. A field star labeled as C (see Figure [TJ) was used 
as a check star, as it was nonvariable and had similar 
brightness to our target. 

Because we did not request observations of standard 
stars for flux calibration in our Gemini program and 
also no standard stars were imaged in the same filter 
by the Gemini South telescope within at least half a 
month before and after our observ ation, we used the BV 
magni tudes of star N measured by I Chevalier fc Ilovaiskvl 
(|1987l ) to obtain absolute magnitudes of the target and 
other stars. The transform ation formula betwee n r' 
and BV magnitudes given bv lFukugita et all (|1996f) was 
used. We found an average r' magnitude of 20.95±0.16 
for 2S 0918—549, where the uncertainty comes from 
the relatively large uncertainties on the BV magni- 
tudes of star N in IChevalier fc Ilovaiskvl (|1987f) . We 
note that with the same transformation, r' = 20.94 for 
2S 0918-549 in IChevalier fc Ilovaiskvl (fl987l ). which in- 
dicates that the binary has not had significant changes in 
its optical brightness. The average r' magnitude of star 
C was 20.98±0.16±0.02, where 0.02 mag is the standard 
deviation of star C measured from 176 frames. 

3. RESULTS 

In Figured] the obtained light curves of 2S 0918-549, 
star TV, and the check star C are shown. A periodic mod- 
ulation in the light curve of 2S 0918—549, while with a 
low-amplitude, is clearly visible. To determine its pe- 
riod, we applied a, phas e-dispersion minimization tech- 
nique (Stellingwcrf 1978) with 16 bins of the full phase in- 
terval (0, 1) used. The resulting periodogram is shown in 
Figure [21 The 9 statistic indicates the detection of a pe- 
riodicity and its two harmonics. Fitt ing the region nea r 
the first minimum with a parabola ( Stclling werJ I1978T) . 
we found period P — 17.4 min. 

In order to quantify the overall periodic modulation in 
the light curve, we fit the light curve with a sinusoid. The 
best-fit has reduced \ 2 = 2-3 for 172 degrees of freedom 
(DOF), and from the best-fit we found P = 17.38 ± 0.13 
min and a semiamplitude of 0.014±0.002 mag. The large 
X 2 value is mainly caused by large scattering of the first 
56 data points due to the poor observing conditions. Ex- 
cluding them and fitting the remaining data points, we 
found reduced x 2 = 1.2 for 117 DOF (P was fixed at 
17.4 min). The obtained semiamplitude was 0.015±0.002 
mag, not having significant changes. Therefore the mod- 
ulation can be described by a sinusoid with a semiampli- 
tude of 0.015 mag. The folded light curve at P = 17 .4 
min as well as the best-fit sinusoid are shown in Fig- 
ure |U The time at the maximum of the sinusoidal fit 
(phase zero) was MJD 54805. 23281±0.00027 (TDB) at 
the solar system barycenter. 

4. DISCUSSION 
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Fig. 2. — r' light curve of 2S 0918—549 (squares). The light curves of the nearby star N (asterisks) and comparison star C (triangles), 
down-shifted by 0.34 and 0.28 mag respectively, are also shown. A sinusoid (dashed curve) is plotted to help indicate the periodic modulation 
detected in the light curve of 2S 0918-549. 



Using Gemini high time-resolution imaging we ob- 
tained an accurate optical light curve of 2S 0918—549 
and have discovered a periodic flux modulation in the 
light curve. A low-amplitude modulation is clearly visi- 
ble and appears to be coherent. Given the known X-ray 
and optical properties of 2S 0918—549, it is very likely 
that we have verified the ultracompact nature of this bi- 
nary and its orbital period is around 17.4 min (see discus- 
sion below). In ad dition to 4U 0513— 40 in the globular 
cluster NGC 1851 dZurek et al.ll2009D and the Galacti c 
disk source 4U 1543-624 (jWang fc Chakrabartvl[200l . 
2S 0918—549 is the third member of the ultracompact 
binaries that have orbital periods around 18 min. 
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Fig. 3. — Phase-dispersion minimization periodogram. The posi- 
tions of the minimum statistic at 17.4 min and its two harmonics 
are indicated by arrows. 

We use the discovered orbital period to estimate the 
mass and radius of the donor. Since the mean den- 
sity of a Roche lobe-filling companion is determined by 



the binary period, our 17 minute period defines a mass- 
radius relation for the companion, shown as the solid 
curve in Figure [5] On the basis of the high Ne/O 
abundance ratio measured through X-ray spectroscopy, 
the donor in 2S 0918—549 has initial ly been suggested 
to be a low-mass C-0 wh ite dwarf (|Juett et all l2001t 
Uuett fe Chakrabartvl 12003) . However, the analysis of 
several type-I X-ray bursts detected from the source 
possibly indi cates that the donor instead is a helium 
white dwarf (|in't Zand et al.ll2005D . In any case to com- 
pare the donor to stellar models, we use the M-R re- 
lations for differen t type s of white dwarfs provided by 
iDelove fc Bildstenl (|2003| ). Because extremely low-mass 
white dwarf donors in ultracompact systems may be ther- 
mally bloat ed compared to cold stars, affecting thei r M- 
R relation (|BildsterJ l200l IDelove fe Bildstenl 12001 . we 
show both cold and hot solutions for pure He, C, and 
O white dwarfs in Figure [5] A helium white dwarf with 
a mass of 0.034-0.039 M and a radius of 0.033-0.035 
R Q , or a C/O white dwarf with a mass of 0.024-0.029 
Mq and a radius of 0.03-0.032 i? Q can fit in the Roche 
lobe-filling donor. 

Modulation of an optical light curve for LMXBs gen- 
erally arises from the companion star that is heated by 
the central X-ray source, with the visible area of the 
heated face varying as a function of orbital phase and 
the superior conjunction of the companion star corre- 
sponding to the observed brightness maximum o f the 
light curve (e.g., Ivan Paradiis fc McClintockl li"995ri . It 
has also been realized that compact LMXBs with ex- 
treme mass ratios (such as ultracompact binarie s) are 
potential superhump sources (jHaswell et al.ll200l . The 
variation in the light curve of a superhump binary arises 
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from an elliptical accretion disk, which is developed when 
the disk extends beyond the 3:1 resonance radius and 
precesses in the iner tial frame due to the tida l force of a 
secondary star (e.g.. lWhitehurst fe Kingl H991). Without 
an i ndependent determination o f the binary period (see, 
e.g. JWang fe Chakrabar tv 2010), we can not distinguish 
between the two possibilities for the modulation seen 
in 2S 0918—549. A super hump modulation may hav e 
an asymmetric shape fe.g. JWang fe Chakrabartvf l2010). 
However in the current light curve we obtained, no asym- 
metry is clearly seen. In either case since superhump pe- 
riods are only a few percent longer than the correspond- 
ing orbital periods, the orbital period of 2S 0918—549 is 
around 17.4 min. 
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Fig. 4.— r' light curve of 2S 0918-549 folded at 17.4 min. Two 
cycles are displayed for clarity. The solid curve indicates the best- 
fit sinusoid with a semiamplitude of 0.015±0.002 mag. The first 56 
data points that were excluded from the fit are shown as crosses. 

Considering that the periodic modulation arises from 
the companion star, the inner face of the compan- 
ion star in 2S 0918—549 is heated by X-ray emission 
from the central neutron star and its effective temper- 
ature can be estimated. The 0.1-200 keV X-ray lu- 
minosity Lx is 1.9xl0 36 d\ erg s _1 , where d§ is the 
source distance assumed to be 5 kpc and the unab- 
sorbed X -ray flux Fx = 6.4 x 10~ 10 erg s" 1 cm~ 2 
given by lin't Zand et~aH (|2005( ) is used. The frac- 
tion / of the X-ray energy absorbed by the compan- 
ion is / = r)*(R 2 /D b ) 2 /4 ~0.004ry,(i? 2 /0.032 P©) 2 [1 + 
(g/0.021)]~ 2 / 3 , where ry» ~ 0.5 is the fraction of the 
received X-ray energy absorbed by the companion, P2 
is the radius of the companion, and Db is the binary 
separation distance (Dt, ~ 1.7 x 10 10 cm for orbital pe- 
riod P or b = 17.4 min). The mass ratio q = M2/M ns , 
where M ns is the neutron star mass and w e assume M2 = 
0.03 M Q and M ns = 1.4 M©. Following lArons fc King] 
(1993), the effective temperature of the companion's in- 
ner face is T = (fLx/TR^a) 1 / 4 ~ 46000dg /2 K, where 
a is the Stefan-Boltzmann constant. The visible area of 
this hot face varies as a function of the orbital phase, 
yielding a modulation of [1 + sini sin(27rf/P rb)] ; where 
i is the inclination angle of the binary (see details in 
lArons fc Kind I1993D . Using such a modulation func- 



tion, we can test how the observed modulation is gener- 
ated. The extinction to the source A y ~ 1.65, estimated 
from A v = A/h/1.79 x 10 21 cm" 2 (jPredehl fc Schmittl 
|1995|) by using hydrogen column density to the source 
N H = 2.95 x 10 21 cm" 2 (|Juett fc Chakrabartv! 120031) . 
By adding a constant flux component (arising from the 
accretion disk) to the modulation function and fitting 
the dereddend light curve of 2S 0918—549 with the mod- 
ulation function, we find i ~ 10°. The estimated low 
inclination angle is consistent with the non-detection of 
orbital signal s at X-ray energie s, which generally indi- 
cates i < 60° ()Frank et al. 2002), and likely explains the 
low-amplitude modulation in the light curve. In order to 
fully explore properties of the binary by fitting the op- 
tical modulatio n, an advanced bin ary light curve model 
is needed (e.g.. iDelove et al1l2008t Wang et al. 2011, in 
preparation) . 
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Fig. 5. — Mass and radius values constrained for the companion 
star in 2S 0918—549. The solid curve is the mass-radius relation 
for a Roche lobe-filling donor in a 17.4 min binary. The dotted, 
dashed, and dash-dot curves represent the model curves f or low - 
mass He, C, and O white dwarfs from Dcloyc & Bildstcn (2003), 
respectively. For each type of white dwarf, both cold (10 4 K) and 
hot (3 xlO 6 K) core temperature model curves are shown. 

As a separate check, we also estimate the distance to 
2S 0918—549. Mass transfer in ultracompact binaries is 
driven by gravitational radiation, and a mass transfer 
rate in 2S 0918—549 can be estimated to be 

M * 6.2 x 1O-M yr-i ( (^-) ' 

/ Porb y 8/3 

\ 17.4 min / 

Since L x = GM ns M/P» = 47rd 2 F x , the unabsorbed 
0.1-200 keV X-ray flux of 2S 0918-549 would imply 
d ^9 kpc, where conservative mass transfer onto a 1.4 
Mq neutron star is assumed and P* ~ 10 km is the 
neutron star radius. The distance value is larger than 
the 4.1-5.4 kpc range derived from type-I X-ray bursts. 
We note that M is sensitive to M 2 - For example, if 
M2 = 0.024 Mq, where an oxygen white dwarf has to 
be assumed, the distance could be lowered to 7 kpc. 
A larger X-ray flux can also help lowered the distance. 
As recorded in an ASCA X-ray observation in 1995, 
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approximately 7 times larg e r X-ray flux was det ected 
(j Juett fc Chakrabartvl [20031 lin't Zand et all [2005h . If 
that ASCA flux is considered, the distance would be 
lowered to ^4 kpc. However, given the known X-ray 
flux history of 2S 0918-549, its 2-10 ke V flux has been 
stable and around 10 -10 erg s _1 cm -2 (jin't Zand et al.1 



2005), suggesting that the ASCA flux was only a one- 
time event. 
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